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ASYMPTOTIC LOCAL EFFICIENCY OF CRAMER VON MISES 
TESTS FOR MULTIVARIATE INDEPENDENCE 
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Deheuvels [J. Multivariate Anal. 11 (1981) 102-113] and Genest 
and Remillard [Test 13 (2004) 335-369] have shown that powerful 
rank tests of multivariate independence can be based on combina- 
tions of asymptotically independent Cramer-von Mises statistics de- 
rived from a Mobius decomposition of the empirical copula process. 
A result on the large-sample behavior of this process under contigu- 
ous sequences of alternatives is used here to give a representation of 
the limiting distribution of such test statistics and to compute their 
relative local asymptotic efficiency. Local power curves and asymp- 
totic relative efficiencies are compared under familiar classes of copula 
alternatives. 

1. Introduction. In a seminal paper concerned with testing the null hy- 
pothesis of independence between the d > 2 components of a multivari- 
ate vector with continuous distribution H and marginals Fi, . . . , Fj, Blum, 
Kiefer and Rosenblatt [1] investigated the use of a Cramer-von Mises statis- 
tic derived from the process 

W n (x) = j #n(z) - J] . x = { Xl ,...,x d )eR d , 

that measures the difference between the empirical distribution function H n 
of H and the product of the marginal empirical distributions Fj^ n asso- 
ciated with the components of the random vector. As Hoeffding [17] had 
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already noted, the asymptotic distribution of this test statistic is generally 
not tractable, and hence tables of critical values are required for its use. 
Such tables were provided by Blum, Kiefer and Rosenblatt [1] themselves in 
the case d = 2. They were later expanded to d > 3 by Csorgo [4] and by Cot- 
terill and Csorgo [2, 3], based on strong approximations of H n . A bootstrap 
approach has also been proposed by Jing and Zhu [18]. 

Despite the anticipation that the Cramer-von Mises statistic fM^dH n 
should be powerful, most subsequent research focused on the case d = 2, 
where alternative tests, typically based on moment characterizations of in- 
dependence, were proposed, for example, by Feuerverger [9], Shih and Louis 
[21], Gieser and Randies [15] and Kallenberg and Ledwina [19]. 

Curiously, the literature seems to have largely ignored a suggestion of 
Blum, Kiefer and Rosenblatt [1] to circumvent the inconvenience caused by 
the complex nature of the limiting distribution of H n . To be specific, let 
X\ = (Xu , . . . , Xid) , ■ ■ ■ , X n = (X n i , . . . , X n( i) be a random sample from the 
distribution H and for arbitrary A C Sd = {1, . . . , d} with \A\ > 1, consider 

1 n 

G A , n (x) = —J^HmXij < Xj ) - F j;n ( Xj )}. 

V n i=l j£A 

Using the multinomial identity 

d 

Y[( a j+bj)=^2 YiajY[bj, ai,...,a d , ,& d eR, 

j=l AcSdjeA jt£A 

with the convention that a product over an empty set of terms equals 1, 
Blum, Kiefer and Rosenblatt [1] showed that H n may be conveniently ex- 
pressed as 

,n\Xj)- 

Ac5 d ,| J 4|>l jeS d \A 

Although their paper focused on the case d = 3, these authors claimed, 
as later confirmed by Dugue [8], that under independence, Ga,u converges 
weakly to a continuous centered Gaussian process with covariance function 

cov A (x,x') = Y[ [Tami{Fj(xj),Fj(x , j )} - F j (x j )F j (x' j )] 
jeA 

and whose eigenvalues, given by ir~ 2 \ A \(ii ■ ■ ■ i\A\)~ 2 f° r *ij---)*'|A| £ ma y 
be deduced from the Karhunen-Loeve decomposition of the Brownian bridge. 
More importantly still, Blum, Kiefer and Rosenblatt [1] and Dugue [8] 
pointed out that the processes Ga,u and Ga',u are asymptotically indepen- 
dent whenever A ^ A' so that suitable combinations of statistics based on 
the individual Ga,u processes could be used to test independence. 
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An obvious limitation of tests based on this approach, however, is the 
dependence of the asymptotic null distribution of the Ga,u on the marginals 
of H . To alleviate this problem, Deheuvels [6] suggested that the original 
observations Xx, . . . ,X n be replaced by their associated rank vectors Ri = 
(Rll, . . . ,Rid), ■ ■■, R n = {Rni, ■ ■ ■ ,Rnd)i where 
n 

Rij = l(Xtj < ±<i<n, l<j<d. 

e=i 

Deheuvels then went on to characterize the asymptotic null behavior of a 
Mobius decomposition of the copula process 

(1.1) C„(«) = Vnjc„(«)-n»ij> 
where the empirical copula [5], defined by 

2 n d 

(1.2) C n (u 1 ,...,u d ) = -YW. 1 ( R ij < n Uj), 

1=1 3=1 

is an estimate of the unique copula C defined implicitly by C{Fi(x\), . . . , 
Fd(xd)} = H(xi, . . . ,Xd) for all x±, . . . ,xa € R. As the latter reduces to C(ui, 
. . . , Ud) = u\ ■ ■ ■ Ud under independence, Deheuvels [6] proposed that this hy- 
pothesis be tested using a Cramer-von Mises statistic based on a decompo- 
sition of C n , but he concentrated his own efforts on the identification of the 
asymptotic null distribution. 

Recently, Genest and Remillard [13] showed how to compute the quantiles 
for the finite-sample and asymptotic null distribution of Cramer-von Mises 
statistics based on C n . Furthermore, they investigated how the 2^ — d — 
1 statistics derived from the rank analogues &A,n of the Ga,u could be 
combined to obtain a global statistic for testing independence. 

This paper enhances the work of Genest and Remillard [13] by comparing 
the power of Cramer-von Mises tests of independence based on the copula 
process C n and of those based on four different combination recipes for the 
terms of its Mobius decomposition. To this end, the local asymptotic behav- 
ior of the copula process C n is characterized in Section 2 under sequences 
of alternatives contiguous to independence. Examples of such sequences are 
considered in Section 3, where it is shown that Clayton's gamma frailty 
model shares with the equicorrelated Gaussian model the surprising prop- 
erty that the limiting behavior of Ga,u is independent of the sequence of 
contiguous alternatives for all A G Sd with \A\ > 2. In Section 4, some test 
statistics for the null hypothesis of independence are considered and their 
asymptotic behavior under contiguous alternatives is specified. Local power 
functions are then computed in Section 5. Finally, Section 6 gives additional 
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comparisons between the test statistics using an extension of Pitman's local 
asymptotic efficiency considered by Genest, Quessy and Remillard [12] in a 
bivariate setting. 

2. Limiting behavior of Cn under contiguous alternatives. Let OcMbe 
a closed interval and C = {Cg : 9 G 0} be a given family of copulas that are 
monotone in 9 with respect to the concordance ordering and for which 9q £Q 
corresponds to independence. Take 5 n — ► S G R such that 9 n = 9q + 5 n / ^Jn G 
G for n sufficiently large. 

Let Q n be the joint distribution function of the random sample (X^ , . . . , 



Xj^ ),..., (^Qi , ■ • ■ , X^j ) from the distribution function Cg n {F\ (x\ ),..., 



Fd(xd)} and let C n be the empirical copula of these observations computed 
using formula (1.2). Finally, let P n be the joint distribution of a sample of 
the same size under the independence distribution F\ x • • • x Fa- 

The limiting behavior of the sequence (C n ) of empirical copula processes 
will be determined under the following conditions: 

(i) Cq is absolutely continuous and its density eg has a square integrable 
right derivative cg(u) = dcg(u)/d9 at 9 = 9q for each u = (u\, . . . ,Ud) G 
(0, l) d , the latter satisfying 



The following result extends Proposition 1 of Genest, Quessy and Remillard 



Proposition 2.1. Suppose that the underlying copula of a given pop- 
ulation belongs to a family C whose members satisfy assumptions (i) and 
(ii) above. Then under Q n , the empirical process C n = \fn{C n — Cg ) con- 
verges in law in T>([0,l] d ) to a continuous Gaussian process C + 8C$ with 
covariance structure A, where for each u,u' G [0,1]^, A(u,u') is given by 




(ii) for every u = (m, . . . , u d ) G (0, l) d , 



one has 




[12]: 



C do (u A v!) + (d- l)Cg Q (u)Cg («') - C 9o (u)C 9o («') ^ ( 




Proof. Introduce L T > nJ = Fj(Xfr } ) and define 
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Also, let Vj, n (u) = £2=i l(E/§° < u)/n. Since B$ = nVj^f), it follows 
from equations (1.1) and (1.2) that 

{d d ~| 

3=1 3=1 > 

Under assumption (i), an application of Theorem 3.10.12 of van der Vaart 
and Wellner [23] implies that under Q n , the sequence (A n ) of processes 
converges in P([0, l] d ) to a continuous Gaussian limit of the form A + 5Cg , 
where Cg is defined as in assumption (ii). 

As a consequence of this result, one has that under Q n , the univariate pro- 
cess A n (l, uj, 1) = ^/n{^ j, n {uj) — Uj} converges in X>([0, 1]) to A(l, Uj, 1) + 
6Co (l,Uj,l) =A(l,Uj,l) since 

C e {l,Uj,l) - C e (l,Uj,l) u-;-Un 

' n n = llm 

+ 6» 0-Qq 0^00 (7 - (7 

Using identities (11) and (12) from Chapter 3 in Shorack and Wellner 
[22], one can see that 

SUp \®j,n{Uj) —Uj\= SUp \^f~ n (Uj) - Uj \ — > 

uje[o,i] tye[o,i] 

in probability for each j G 5^, whence y/n{^J n (uj) — Uj} —>■ —A(l,uj, 1) in 
T>([0, 1]). Finally, the second summand in (2.1) can be rewritten as 

EfnV){ n *,>(«/)}^w>*)-«*}. 
fe=i\i=i / ij=fc+i j 

Accordingly, one may conclude that under Q n , the process C n converges in 
D([0, l] d ) to C + <5C\,, where 

d 

C(u) = A(u) - J2 M 1 , u k, 1) n U 3 

k=l j^k 

is the limiting process of C n under P n that was identified by Ganssler and 
Stute [10]. Finally, a straightforward computation shows that the limiting 
covariance function of C n is A, given that 5Cg is a deterministic term. □ 

Deheuvels [6] proposed the decomposition of <C n into a collection of asymp- 
totically independent, centered Gaussian processes having a simple covari- 
ance function under the null hypothesis of independence. Specifically, for 
each A C S4, define the linear operator Ma such that 

1 " 

Ga,„(«) = M A {C n (u)} = -= II i^ R v ^ nu i) ~ u i) 

v n i=i je A 

= J2 { -l)\A\B\ Cn{u B )Ceo{u A\B^ 
BcA 



C eo (l, Uj ,l) = lim ^ ^ ^ = lim ^ = 0. 
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where in general, u B = (uf , ■ ■ ■ ,u d ) with equal to Uj when j G B and 
equal to 1 otherwise. The result is that under P n , one obtains a decomposi- 
tion 

C n (u) = G A,n( U ) II U 3 

AcS d ,\A\>l jeS d \A 

of the empirical process C n into 2 d — d — 1 subprocesses that converge jointly 
to a vector of continuous, centered Gaussian processes Ga = Ma(G). Fur- 
thermore, their asymptotic covariance structure is given by 

T AA ,(u,u') = cov{G A (u),G A ,(u')} = 1(A = A')H Jiu^u'j), 

j£A 

where 7(u,u') = u A u' — uu' . In other words, the G^ jn are asymptotically 
independent and their limiting covariance is the same as that of a product 
of independent Brownian bridges. 

The next result, which is a straightforward consequence of Proposition 2.1, 
gives the asymptotic representation of Ga,u under Q n . 

Corollary 2.2. Let C be a given family of copulas whose members 
satisfy assumptions (i) and (ii) above. Then under Q n , the empirical pro- 
cesses Ga,u converge jointly in T>([0, l] d ) to a vector of continuous Gaussian 
processes Ga + $Pa with covariance structure Ta.A', where fJ, A = -M-a(C9 ). 

3. Examples. The drift term /j,a(u) identified in Corollary 2.2 can be 
computed explicitly in many families of copulas. The following result, known 
as the Mobius inversion formula, will be useful for that purpose. 



Lemma 3.1. Let f be a function defined on the subsets B of S d . For any 
AcS d , setF(A) = J2BcAf(B). Then f{A) = Ebca(-^ B{ F(B). 

3.1. The multivariate equicorrelated normal copula. Let $ denote the 
cumulative distribution function of a standard normal random variable. The 
multivariate normal copula with d x d correlation matrix S = ((Tjk) is defined 
by C s (u) = H Tl {$- 1 (u 1 ),...,<5>- 1 (u d )}, where 



■>:\ 



1 



exp(-y T S 1 y/2)dy d ---dy 1 



, (27r) d / 2 |S| 1 /2 

Consider the equicorrelated case in which o~jj = 1 for all 1 < j < d and 
a ik = P f° r au J k- Write H p = Hy, and H p = dH p /dp. Then 



dp 



and 



p=0 



ld_ 

2 dp* 



p=0 j<k 
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where the latter identity follows from the fact that 



T V -1 0-~P) 
x h x = 



d-2 



Now if <p(t) = d®(t)/dt, one gets 



d 

{{d - 2)p + 1} x i ~ 2 P E x i x k 

1=1 j<k 



lim H p (x) = j ■■■ f ( ^ yjy k ) J J| </?(y^) I dy d ■ ■ ■ dyi 

=e{ n ^l^M^i^oWE^^^ 



Thus, 



, <£(x,) $(xk) 



C 9o (u) = lim V), • • • , <S>-\u d )} 



An application of the Mobius inversion formula then yields 

^(n)=i(i^i=2)n^" i (^)}- 

Hence, when looking at the Mobius decomposition of the multivariate 
normal model, only the Ga,u with \A\ = 2 have a limiting distribution that 
differs under the null hypothesis and contiguous alternatives. Accordingly, 
tests of independence should be based only on the latter because the in- 
clusion of functions of Ga,u for any \A\ > 2 would contribute nothing to 
the overall power of the procedure. This observation does not come as a 
total surprise, given that the multivariate Gaussian dependence structure is 
completely characterized by the pairwise interactions among the variables. 

3.2. One-parameter multivariate Farlie-Gumbel-Morgenstern copula. Ad- 
variate version of this system of copulas is defined for 6 £ [—1, 1] by 

d 

C e (u) = C eo (u) + eY[ Uj (l - Uj ), ue [0,l] d , 
i=i 

with #o = corresponding to independence. It follows easily that Cg (u) = 
u\{\ — u\) x • • • x Ud(l — Ud). Since Cg (u B ) = whenever B j^Sd, one gets 

fi A (u) = (-l) d C eo (u)l(A = Sd). 

Thus, in this case, [ia vanishes unless A = Sd, implying that in contrast to 
the multivariate normal case, tests of independence based on the G A ^ n with 
any | vl| < d would have no power in the neighborhood of independence. 



8 C. GENEST, J.-F. QUESSY AND B. REMILLARD 

3.3. Archimedean copulas. Following Genest and MacKay [11] and Nel- 
son [20], a copula is called Archimedean whenever it can be expressed in the 

form C(u) = <j)~ l {4>{u\) H \-4>(ud)}, for some generator (ft : (0, 1] — > [0, oo). 

To insure that C is a copula, it suffices that 4>{1) = and (—iyd^(j)~ 1 (t)/dP > 
for every j 6 Sd- 

Proposition 3.2. LetC be a parametric family of Archimedean copulas 
with generator <j)Q such that 4>e(t) — logi and 4>'a(t) — ► — 1/t as 9 — > 9q. 
Further, assume that C satisfies assumptions (i) and (ii) stated in Section 2. 
Then 

(3-1) 7^777] = K{Ce (u)} - X>oK0 

and 

(3.2) n A (u) =C dQ {u A ) ^ (-l)^<j>e {C 9o (u B )}. 

BcA 

Proof. Noting that (fio{Cg(u)} = 4>g(ui) + • • • + 4>$(u ( i) and applying the 

chain rule, one finds 4>e{Ce(u)} + Ce(u)(j)' d {Cg(u)} = <pe(ui) H h 4>e(ud). 

Equation (3.1) follows by taking the limit as 9 — > #o- As for (3.2), it follows 
from substitution of Cg into the formula for Ma, combined with the fact 
that when \A\ > 2, the Mobius inversion formula yields 

e E(- i ) |AB| ^oK')=i(i^i=i)n^o(%)=o. □ 

BcAjeB j£A 

Example 3.1. Assumptions (i) and (ii) can easily be checked for Frank's 
family of ti-variate copulas generated by 4>e{t) = \og{{e~ e — l)/{e~ et — 1)} for 
t £ (0, 1] with 9 £ [£ d , oo), where -oo = £ 2 < £3 < ■ ■ ■ < 4o = 0. Here, 9 = 
corresponds to independence and (f>6 {t) = (t — l)/2. Proposition 3.2 and the 
Mobius inversion formula together yield 

Ce (u) = hj 6o (u) jd - 1 + C do (u) - E Uj j 

and 

M«) = \c d0 {u A ) E (-^'{(^(u*) - 1} 



= kfcCti^) e (-i)^ b| q (o = i n - !)■ 

Z BCA Z j€A 

Here, /U^ 7^ for every AcSj with \A\ > 2. 
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Remark 3.1. The <i-variate Ali-Mikhail-Haq Archimedean system of 
copulas with 9 £ [0, 1] is generated by (j> (t) = (1 - 9)^ \og{9 + (1 - 9)/t} for 
t £ (0, 1] and 9q = Q corresponds to independence. For this system, one gets 
(j)0 Q {t) =t—l — logt and hence ha{u) = Ylj£A u j( u j ~ 1)) which is the same 
as for Frank's family, up to a multiplicative constant. A similar conclusion 
holds for another version of the Farlie-Gumbel-Morgenstern copula, namely 

C e {u) = C do (u) + 9C 9o (u) j d - 1 + C 6a («) - J2 u 5 j . 



Example 3.2. Assumptions (i) and (ii) can also be easily verified for 
Clayton's (i-variate family of copulas whose generator, defined for all 9 £ 
[0, oo), is given by 4>e(t) = (t~ e — l)/9 for t £ (0, 1]. Note that independence 
corresponds to the value 9q = 0. One easily finds that (j)e {t) = (logi) 2 /2. 
Calling on Proposition 3.2 and the Mobius inversion formula, one obtains 

Co (u) = Cg (u) log u k 

j<k 

and 

(J,a(u) = H\A\ = 2) Y[ uj logitj. 



Remark 3.2. The dependence structure induced by Clayton's copula 
is also known as the gamma frailty model. It may come somewhat as a 
surprise that in this case, Ha(u) = unless \ A\ = 2. In other words, Clayton's 
copula shares with the multivariate normal model the property that tests of 
independence based on terms n of the Mobius decomposition with \A\ > 2 
would have no power whatsoever in the neighborhood of independence. 

The Gumbel-Barnett system of copulas provides another example of this 
curious phenomenon. Copulas in this class are generated by (f>g(t) = log(l — 
9 logt)/ 9 for t G (0, 1] with 9 £ [0, 1] and 9$ = corresponding to indepen- 
dence. A simple calculation shows that 4>e (t) = — (logt) 2 /2, and hence the 
formulas are the same as for Clayton's copula, up to a change in sign. 



Example 3.3. Assumptions (i) and (ii) can also easily be verified for 
the so-called Gumbel-Hougaard family of copulas whose generator is defined 
for 9 £ [0,1) by cp e (t) = |logt| 1/(1 ~ 9) for t £ (0,1]. In that case, one finds 
4>g (t) = — (logt) log (log 1/i). Thus, Proposition 3.2 entails that 

{d / d 

- log uj log Y lo S u k/ lo g u i 
i=x \k=i 
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and 

H A (u) = -C 9o (u A ) J2 ( E lo g^) lQ g(- E lo ? 

BcA ^jeB ' ^ j£B 

for all u G (0, l) d . Here, again, ha^O for every ic5d with |A| > 2. 

4. Limiting distributions of Cramer von Mises functionals. In the ab- 
sence of information about the marginal distributions of a multivariate pop- 
ulation, a valid testing procedure for independence should be based on some 
version of the empirical copula process C n . To improve convergence and 
reduce bias in finite samples, a centered version of C n will be used in the 
sequel. The latter is defined by 

C n ( U ) = V^\c n {u) ~ J] Cn(«i)|, 

where U n is the cumulative distribution function of a uniformly distributed 
random variable on the set {1/n, . . . , n/n}. It is clear that C n and C n have 
the same limiting behavior, so the asymptotic results of Section 2 also apply 
to C n and henceforth to Ga,u = A^A(C n ) for Ac Sd- 

A natural way to test for independence is to consider a global measure 
of discrepancy computed from C n or from a combination of distances com- 
puted for each of the &A,n taken individually. Obvious candidates are the 
Kolmogorov-Smirnov statistics 

S n = sup |C n (u)| and Sa,u = sup \Ga,ti(u)\ 

uG[0,l] d ue[0,l] d 

and the Cramer-von Mises functionals 

B n = {C n (u)} 2 du and B A , n = {GA,n(u)} 2 du. 
J(0,l) d J{0,l) d 

In the sequel, however, attention will be limited to B n and Ba,u- Note 
that these statistics can be expressed as functions of the ranks through 

1 * x-^ 1-r / . rCik V Kjk 



i=i j=i fe=i 

'n(n- 1) -R ik {R ik - 1) 



i=ifc=i 



+ n 



(n-l)(2n-l) 
6n 2 



2ra 2 



d 
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and 



where 



Y n n 

BA,n = — X] X] II B> n {Rik,Rjk), 
i=l j=l keA 



(n+l)(2n + l) s(s-l) t(t - 1) sVt 
D n(s,t) = — -„ + „_o + 



6ra 2 2ra 2 2n 2 



4.1. Asymptotics for B n . Since no explicit Karhunen-Loeve expansion 
for C is available when d > 2, even under P n , the asymptotic null distribution 
of B n cannot be computed analytically. This is due to the unwieldy form of 
the covariance function A. 

If C is a family of copulas whose members satisfy assumptions (i) and (ii), 
one has from Proposition 2.1 that B n converges in law under Q n to 



(4.1) B= / {C(u) + 5C eo (u)Y du, 

J(o,i) d 

which is a Cramer-von Mises functional of a Gaussian process with mean 
5Ce Q and covariance function A. In order to approximate the distribution 
of B, a procedure due to Deheuvels and Martynov [7] is adopted here. Specif- 
ically, the proposed approximation is 

-. m 

(4-2) B = -£4 2 , 

m 7~r 

where £ = A({7) + V{U)Z is an m-variate vector constructed from two in- 
dependent vectors U and Z whose components are independent and U(0, 1) 
and JV(0,1), respectively. Here, A(m, u m ) = S(Cg (m), C 6o (u m )) T 
and V(u\, . . . ,u m ) is the Cholesky decomposition of the covariance matrix 
£ with components E,-£ = A(v,j,Uk) for j, k s5 m . Deheuvels and Martynov 
show that the variance of the approximation error is 0(l/m). 



4.2. Asymptotics for Ba,u- In view of the fact that the covariance struc- 
ture of Ga is a product of covariance functions associated with Brownian 
bridges, it follows from standard theory [22, p. 213] that under P n , the 
limiting process Ga admits the representation 

(4.3) G A (u)= Yj \J~^~t z ifi( u )i n / = (lj)jeA, 

where the Z 7 are independent AT(0, 1) random variables and 

A 7 = II (^i)" 2 ' AC") = II V^sinfrj-Traj), u £ [0, l] d . 
j&A jeA 
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An idea of Deheuvels [6], later exploited by Genest and Remillard [13], is 
to base a test of independence on some combination of the asymptotically in- 
dependent statistics Ba u - In view of Corollary 2.2, the limiting distribution 
of Ba,u under Q n is given by 

{Ga(u) + $Ha{u)} 2 du. 

l(0,l) d 

As already shown by Deheuvels [6], formula (4.3) implies that 
[ {G A (u)} 2 du= Vv ^7=IlV 

It follows that 

M a = Qa + 2S J2 A 7 / 7 ,a^ 7 + 5 2 /4, 

7 eNl- 4 l 

where 

I A = [ {n A (u)} 2 du and I-., a = —~= [ fiA(u)fj(u)du. 
J{o,i) d \/\ J(o,i) d 

The limiting distribution of Ba,u under Q n is given below. The result 
follows from direct substitution of Parseval's identity I a = Z) 7 gnI a I A7Z 2 , A 
into the integral representation of . 

Proposition 4.1. If C is a family of copulas whose members satisfy 
assumptions (i) and (ii), then the asymptotic distribution of Ba,ti under Q n 
is given by the weighted sum 

m a= E A 7 (Z 7 + 5/ 7 , A ) 2 . 
7 gnI a i 

Note that there also exists a representation like (4.3) for C, but its weights 
are unknown. Nevertheless, their sum is finite. Once again calling on Parse- 
val's identity and using (4.1), one obtains 

B= J2 A 7 (^ 7 + <5/ 7 ) 2 , 

7SN d 



where 



and 



7 A 7 = / / T(u,u') du' du < oo 

^N" J (0,1)1 J (0,1)1 



V A 7 / 7 = / {Ce (u)} 2 du<& 
J{o,i) d 
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The unknown quantities A 7 and 7~ could be approximated numerically. 

For many systems of distributions, the drift term has the simple form 
Ha(u) = Y\j & At JL { u j)- I n that case, it is easy to see that 

2 \A\/2 ,1 

(4.4) Iry A = — ] [ f(jj), where f(k)= / n{u) s\n(kiru) du. 

Example 4.1 (Equicorrelated Gaussian copulas). For this model, (4.4) 
and results from Section 3.1 together imply that I 7> a = 1(1^1 = 2)27r 2 IljeA lj9{lj)i 
where g{k) = J ^j^" 1 (u)} sin(kTru) du. 

Example 4.2 (Farlie-Gumbel-Morgenstern copulas). For this model, 
results from Section 3.2 imply that I 7|j4 = {-l) d l{A = S d )2 M / 2 A 7 if 71, . . ., 
jd are all odd and I yj A = otherwise. 

Example 4.3 (Frank and Ali-Mikhail-Haq copulas). For these mod- 
els, calculations based on material from Example 3.1 lead to I 7jJ 4 = (— l)'" 4 ' 
2 5 l j4 l/ 2 ~ 1 A 7 if 7j is odd for every j £ A and I~ a = otherwise. 

Example 4.4 (Clayton and Gumbel-Barnett distributions). For these 
models, the observations already made in Example 3.3 and Remark 3.2 yield 
I 7jA = 1(\A\ = 2)2^- 2 Y\ j&A SI{ lj Ti)h j , where SI(x) = /* t" 1 sin(i) dt. 

4.3. Combination of independent statistics. A simple test of indepen- 
dence consists of rejecting the null hypothesis whenever the observed value 
of B n exceeds the (1 — a)-percentile of its asymptotic distribution under 
P n . The critical value ^(a) such that P{B n > qB{ot)\P n } — ► ct as n — > 00 
is easily approximated using formula (4.2) with 5 = 0. For a = 0.05, these 
values are to be found in Table 1. 

However, potentially more efficient methods could be based on combi- 
nations of the asymptotically independent -Ezra's with A £ Sd- Four such 
procedures are considered here. The first is inspired by Ghoudi, Kulperger 
and Remillard [14], while the second and third were studied by Genest and 
Remillard [13]. 

4.3.1. Linear combination rule. Base the test on 

L n = Ba,h- 
\A\>1 

Table 1 gives approximate values for (?z,(a) such that P{L n > qL(ot)\P n } — ► ct 
as n — > 00. 
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Table 1 

Approximate critical values of the tests based on B n , L n , W, 
M n and T n , at the 5% level 





d = 3 


d = 4 


d = 5 


g B (0.05) 


0.05669 


0.04124 


0.02549 


gi,(0.05) 


0.14045 


0.26002 


0.42046 


g^(0.05) 


18.50403 


50.54507 


134.38756 


ga(a') 


0.08518 


0.10126 


0.12731 


93 (a') 


0.01006 


0.01186 


0.01326 


94(a') 




0.00159 


0.00175 


95 (a') 






0.00022 


g T (0.05) 


15.31231 


35.09049 


71.00888 



4.3.2. Dependogram method. Base the test on 

M n = max.{B A ,n/q\A\(oi')}, 
\A\>1 

where q\A\(ct) is such that P{Ba,u > Q\A\( a )\Pn} — ► ct as n — > oo and where 
(4.5) a ' = i_(i_ a )i/(2 d -^i) 

is chosen so that by performing each test at level a' , the global level of the 
procedure is a. The first few simulated critical values q\M are reported in 
Table 1. Note that P(M„ > l\P n ) -> a as n -> oo. 

4.3.3. Fisher's approach. Base the test on 

T n = -2^ log{l -i^n^n)}, 

|A|>1 

where F^ n {x) = P(B J 4 ) „ < x). Under P n , T n converges in law to 

-2 £ log{l-<^(Q A )}, 

|A|>1 

which is chi-square with 2{2 d -d-l) de grees of freedom. Hence, the critical 
value of a test based on T n is given by 

q T (a) = K _1 (1 - a), where K(x) = P{X2(2 d -d-i) - x i- 

Table 1 reports the asymptotic values of q"r(0.05) for d E {3,4,5}, as per 
Genest and Remillard [13]. 
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4.3.4. Weighted linear combination rule. Base the test on 

W n = E ^ B A.n. 

\A\>\ 

Table 1 gives approximate values for qw( a ) such that P{Vt / „ > qw(a)\P n } — > 
a as n — > oo. This procedure is inspired by Proposition A.l since in view 
of the latter, — 21og{l — FA, n (BA,n)} ~ ^ 2 ^Ba,u whenever Ba,u is large. 
Therefore, W n should approximate T n well under fixed alternatives Cq be- 
cause in that case, BA, n / n tends to a positive constant. 

5. Comparison of local power functions. One way to compare compet- 
ing test procedures for independence is to evaluate their asymptotic power 
function in a neighborhood of 9 = 9q, that is, under copula alternatives C$ n 
that form a contiguous sequence (Q n )- Specifically, let S n be some statistic 
for independence with asymptotic critical value qs(ct). The associated local 
power function is defined as 

P 8 (a,S)= lim P{S n >q s (a)\Q n }. 

Analytic expressions for this function are given in Section 5.1 for the 
combined test statistics L n , M n and W n . No similar form could be obtained 
for B n or for Fisher's procedure T n ; however, see Remark 5.1. For numerical 
comparisons, see Sections 5.2 and 5.3. 

5.1. Analytic expressions for local power functions. A result that will 
prove useful in the sequel is the formula of Gil-Pelaez [16], which says that 
if X is a random variable with characteristic function (ftx, then 

1 1 r°° 

P(X>x) = - + — lm{r l e- ixt <t>x{t)}dt, 

I ITT J-oo 

where Im(^) denotes the imaginary part of any complex number z. 
To use this identity in the present context, let 

^^) = (l-2zt)- 1/2 ex P ( T ^) 

_ _|_ ^2^1/4 e -2t 2 <5 2 /(l+4t 2 ) e iarctan(2t)/2+it5 2 /(l+4t 2 ) 

be the characteristic function of a noncentral chi-square variable (Z + 5) 2 , 
where Z~AA(0, 1). From Proposition 4.1, it follows that 

0B A (M) = E(e^)= J] ftiX^SI^A) 

_ c U\ -25 2 t 2 n A ,i(t) jK A ,2(t)+i8 2 K A , 3 (t) 
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where 



U(t)= II (l + 4t 2 A?)- 1 /4 ) KA>1 (t)= £ A2/2 A /(i + 4 t 2 A2 



«A,2W = ^ £ arctan(2U 7 ), K Al3 (t)=t £ A 7 / 2 A /(l + 4i 2 A 2 ). 

7 eNl^l 7 gNI a I 

According to Proposition 4.1, the limiting distribution of L n under Q n is 
L = J2\a\>i^a, so its asymptotic characteristic function is given by 

\A\>1 

where £(t) = II|A|>i£a(*) and = E|A|>1 «A,i(£) for i E {1,2,3}. 

An application of the Gil-Pelaez formula then yields the following result: 

Proposition 5.1. If C is a family of copulas whose members satisfy 
assumptions (i) and (ii) ; then under Q n , /3^(a,<5) = P{L > q^^a)} , where 

1 1 Z" 00 sin{«(x,t)} , 
s dt 



1 1 f 

P(L>x) = - + - / 

2 7T Jo 



tat) 



with 

|A|>1 7G N 

and 

5 2 t 2 _ _ A, 2 / 2 



«0M) = -y + ^ E E jarctan(A 7 t) + 5 2 ]_ ^ | 



c(t)=«p(*f e E n n a+^ 2 



2^1/4 

A|>1 7 gNI a ! " ' ' >7 " 7 |A|>l 7 eNl A l 



A similar result for the test based on M n follows from the fact that 

Jim P(M n >l) = l- [] P{BA<g|A|(«')} 
|A|>1 

(5.1) =1- I] {1-Pa(o/,S)} 

\A\>1 

with a' defined as in (4.5). 



Proposition 5.2. Let M be the limit in distribution of M n under Q n 
and define (3m to be its associated local power function. If C is a family 
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of copulas whose members satisfy assumptions (i) and (ii), then under Q n , 
(3m(c<,8) is given by (5.1) and for all x > 0, 



2 Wo *Ca(*) 
ka(x,*) = — 2^ jarctan(A 7 t) + (5 — ^2 f 

and 

a(t)=exp(^ E n (i+W 4 - 

Next, one can see that the limiting distribution of W n under Q n is W = 
2~Z|A|>i 7r2 '^'®A, so its asymptotic characteristic function is given by 

|A|>1 

where for z G {1, 2, 3}, 

m*)= n ^(^ 2|A| ) and «i.w(<)= e ^(^ 2|A| )- 

|A|>1 |A|>1 
The following proposition gives the local power function of W: 

Proposition 5.3. If C is a family of copulas whose members satisfy 
assumptions (i) and (ii), then under Q n , (3w{ol,8) = Y*{W > qw(o)}, where 

1 , 1 f°° am{K W (x,t)} 



1 If 

P(W>x) = - + - \ 

2 IT Jo 



Kw(t) 



with 

K W (x,t) = -y + ^ E E (arctan^Tr^lt) + 5 2 ^xl ' W } 

and 

/J?2+2 \2„4|A|r2 N 

cm=«p(*f e e ^ n n d+w. 

v " |A|>l 7 gNl A l T |A|>l 7 eNl^l 
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\A\>1 



Remark 5.1. For Fisher's test, no explicit representation for /3r(ct,5) 
seems possible. Under Q n , however, T n converges in distribution to 

-2 £ log{l-gj^(B A )} 

where 

k W>-2~*J ~^CoW 

with 

xt 1 

K (x,t) = -- + - arctan(A 7 t) and = JJ (1 + *?* 2 ) 7 - 

7< =Nl A l 

An approximation is then given by 

1 N 

P T (a,5) = -Y / HT i >q T (a)}, 

where 



iV . 
i=l 



T i = -2J2 log{l-^|(JEM} 

\A\>1 

and for iV sufficiently large, i, . . . , Man are mutually independent obser- 
vations from a finite-sum approximation of B^. 



5.2. Power comparisons for statistics involving all A G Sd with \A\ > 1. 
Figure 1 graphically compares the power of five different tests for trivariate 
independence based on statistics B n , L n , M n , T n and W n . These compar- 
isons were carried out at the 5% level for four parametric classes of three- 
dimensional copula alternatives considered in Section 3, namely (a) the 
equicorrelated Gaussian copula, (b) the Farlie-Gumbel-Morgensten copula, 
(c) the Frank or Ali-Mikhail-Haq copula and (d) the Clayton or Gumbel- 
Barnett copula. In view of the considerations made in Sections 3 and 4, 
the choice of representative within each class is irrelevant for the following 
asymptotic local power comparisons. 

The local power curves corresponding to L n , M n and W n were evaluated 
by means of numerical integration of the formulas given in Propositions 5.1, 
5.2 and 5.3. In each case, the integral was computed by the trapezoidal rule 
on a domain of the form [0, K] for suitably large K; all infinite sums were 
truncated to 40 terms in each index. This was sufficient to ensure numerical 
accuracy. As for the local power curves associated with B n and T n , they 
were obtained via Monte Carlo simulation, using 10,000 repetitions for each 
point on the curve. 
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deiia delta 

Fig. 1. Comparative local power of B n ( — ), T n ( — ), M n (■■■), L n ( — ) and W n 

( ) for four different classes of trivariate copula alternatives: the equicorrelated Gaussian 

(upper left panel), the Farlie-Gumbel-Morgenstern (upper right), Frank/ Ali-Mikhail-Haq 
(lower left) and the Clayton/Gumbel-Barnett copula (lower right). 



Looking at Figure 1, one can see that the test based on B n is best for 
the Clayton/Gumbel-Barnett models in the neighborhood of independence. 
It is also close to optimal for the Frank/ Ali-Mikhail-Haq models, but its 
behavior is much less satisfactory for the Farlie-Gumbel-Morgenstern copula 
and especially for the equicorrelated Gaussian structure. 

Putting B n aside, the smallest local power is yielded by M n for all struc- 
tures considered except the Farlie-Gumbel-Morgenstern. In the latter model, 
L n is the least powerful locally. This results from the fact that this procedure 
equally weights the four Ba,u statistics of size \A\ = 2 or 3, whereas only the 
latter has any power in detecting dependence, as observed in Section 3.2. 
In contrast, note the excellent performance of W n , which weights the Ba,u 
proportionally to -/r 2 '" 4 '. 

Given that W n is an approximation of T n , as seen in Section 4.3.4, it is not 
surprising that they behave similarly in all cases. While Fisher's procedure 
T n ranks first in the trivariate normal model, there does not seem to be 
much difference between T n and W n in the other models considered. 
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Table 2 

Approximate critical values of the tests based on in, 2, M n ^ 
and T n ,2) at the 5% level 





d = 3 


d = 4 


d = 5 


<7l 2 (0.05) 


0.13562 


0.23917 


0.36848 


ga(a") 


0.07479 


0.09020 


0.09714 


?t 2 (0.05) 


12.20343 


20.32429 


29.63573 



5.3. Power comparisons for statistics based on sets A with \A\ = 2. In 
Section 3, it was seen that for Clayton/Gumbel-Barnett and equicorrelated 
Gaussian contiguous alternatives, ha = when \A\ > 2. As a consequence, 
the asymptotic distribution of Ba.u under Q n is then the same as that 
under P n . A loss in efficiency may thus be expected to occur when testing 
for independence within these models whenever a test statistic combines all 
possible E>A, n , rather than only those for which |A| = 2. 

For such alternatives, possibly more efficient tests could be based on 

L n ,2= B A , n , M n 2 = maxi ^ A '" 
|AHJ \A\=2{q 2 (a") 

and 

T n ,2= \og{l-F A> n{B A> n)}, 

\A\=2 

where a" = 1 — (1 — a) 2 ^ d<yd ~ 1 ^ and 92 (a") is given in Table 2. Analytic 
expressions for the asymptotic local power curves of these statistics can 
be derived from straightforward adaptations of Propositions 5.1, 5.2 and 
5.3. Table 2 gives the critical values, at the 5% level, of the tests based on 
these statistics, based on Monte Carlo simulations. These values were used 
in Figure 2 to compare the local power curves of these statistics to their 
analogous versions based on all \A\ > 1. 

In the case al = 3, only the statistic Ba,u with |^4| = 3 has no local power in 
testing for independence in the equicorrelated Gaussian model. The loss of 
power caused by its inclusion is most apparent for the dependogram statistic 
M n , as illustrated in the left panel of Figure 2. The right panel of the same 
figure shows that the loss is much less dramatic (if indeed there is any) for 
T n compared to T n> 2- This was perhaps to be expected, considering that 
Fisher's test based on T n is probably close to optimal in this context. 

6. Asymptotic efficiency results. Given two competing test statistics 
S nt i and S n> 2, let 

/3si(a,6)= lim P{S n ^ > q St {a)\Q n } 
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FlG. 2. Comparison of the local power of two test statistics for independence in the tri- 
variate, equicorrelated Gaussian model. Left panel: M n> 2 (solid line) versus M n (broken 
line). Right panel: T n> 2 (solid line) versus T n (broken line). 



and suppose that {(3s t {ot,5) — a}/5 — > as S — ► for i G {1,2}. Genest, 
Quessy and Remillard [12] argue that an appropriate measure of the asymp- 
totic relative efficiency of S n< i with respect to 5 nj 2 is then given by 

f( ..s r \ v Psi(a>,5)-a \im s ^ Q {f3 Sl (a,5) - a}/5 2 
(6.1) e 12 {a) = hm — - - — = - — — — — 

5^o (3s 2 {a,5) - a hm 5 ^ {/3s , 2 (Q,5) - a}/5 2 

which is a limiting ratio of curvatures of local power functions at 5 = 0. 
Analytic expressions for these curvatures can be obtained for the tests of 
independence based on L n , W n and M n . They are given in the next propo- 
sition, whose proof depends on Proposition A. 2. In the sequel, 0a denotes 
the local power function associated with the statistic Bj\ n . 

Proposition 6.1. Let h w \M, h Wj i and h W) w be the densities of 

Qa + wxI J2®A + wxl and J2 ^ AI Qa + w X l 

\A\>1 \A\>1 

respectively, where in each case, the summands are taken to be independent. 
Then as 5^0, one has 



(ii) 



(iii) 



{P A (a,5)-f3 A (a,0)}/6 2 ^ £ ^^^{q^a)}; 

{(3 L (a,5)-f3 L (a,0)}/5 2 ^ E ^,aK,lUl(ol)}; 

|A|>i 7 eNl j4 l 

{(3 w (a,S)-(3 w (a,0)}/6 2 ^ J2 E V 2|A|/ 7, ^A 7 ,w{<7w(a)}; 

|A|>i 7 eNl A l 
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(iv) given a' defined as in (4.5), 

{/3 M (a,S) - l3 M (a,0)}/5 2 -» ( — ^ ^ A 7 I^ A /i A7i | A |{g| A |(a')}- 

Proof. Recalling that 

IU= E A 7 (Z 7 + 5/ 7iA ) 2 , L= ^ B A and TV = ]T tt 2 ^M a , 
7 GNl A l |A|>1 |A|>1 

one can derive results (i), (ii) and (iii) easily using Proposition A. 2. Now, 
using (i), one has /3a(oc,5) = Pa(&,0) + 5 2 £a(&) + o(<5 2 ), where 

7 eNi A l 

Statement (iv) then follows, on writing 

(3 M (a,6)-f3 M (a,0)= JJ {1 - (3 A (a',0)} - ]J {1 - p A (a ,5)} 

\A\>1 \A\>1 

= «5 2 ]T£ A (a') II {l-PD(a',0)} + o(5 2 ) 

\A\>1 \D\>1,D^A 

5 2 U(a')+o(5 2 ). □ 

|A|>1 

The above result, combined with equation (6.1), makes it possible to de- 
termine the local asymptotic efficiencies of statistics L n , M n and W n . These 
values are given in Table 3, along with those corresponding to L n ^ and M n ^, 
for which a simple adaptation of Proposition 6.1 must be used. 

The results in Table 3 are generally in line with those reported in Sec- 
tion 5, bearing in mind that due to the lack of analytical expressions for 
their local power curves, the Cramer-von Mises and Fisher test procedures, 
based respectively on B n and T n , could not be included in the efficiency 
comparisons. 

The following observations are offered as concluding remarks. 

(a) Among trivariate tests of independence based on L n , M n and W n 
which involve all of the Ba,u , L n was most efficient and M n was least efficient 
in all models considered except the Farlie-Gumbel-Morgenstern class of 
copulas. In the latter case, M n was by far the best choice. 

(b) Tests based on L n ^ or M n ^ were totally inefficient when dependence 
entered through Farlie-Gumbel-Morgenstern trivariate copulas. This was 
to be expected because in this case, only the term involving Ba,u with 
\A\ = 3 has a limiting distribution that differs under the null and under the 
alternative. Tests based on L n and W n hardly did any better. 
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Table 3 

Local asymptotic relative efficiency of Cramer-von Mises tests for trivariate independence 



Alternative Best statistic among Relative efficiency of the statistic 



trivariate 


L n , M n , W n , L n ,2, M„, 2 




L n ,2 




M„, 2 




model 


for the chosen model 




with respect to the best 




Equicorrelated 


L n ,2 


98.56 


100.0 


44.66 


88.78 


79.45 


Gaussian 














Farlie- 


M n 


3.71 





100.0 





32.95 


Gumbel- 














Morgenstern 














Frank Ali- 


L n ,2 


99.34 


100.0 


65.28 


88.92 


86.09 


Mikhail-Haq 














Clayton 


L n ,2 


98.55 


100.0 


43.27 


87.21 


79.79 


Gumbel- 














Barnett 















(c) For the other models, a loss in efficiency was observed when going 
from L Ut 2 to L n and from M n ^ to M n . In the equicorrelated Gaussian and 
Clayton/Gumbel-Barnett cases, this was expected since it was argued earlier 
that the inclusion of statistics Ba, u with \A\ > 2 is then likely to dilute the 
power of the overall procedure. An explanation for the same occurrence in 
the Frank and Ali-Mikhail-Haq models is still lacking. 

If nothing else, this study provides a new illustration of the truism that 
no single procedure could ever be declared best for simultaneously testing 
for all forms of multivariate dependence. More importantly, however, the 
results reported herein shed new light on unsuspected commonalities among 
classes of dependence models that may be superficially perceived as quite 
different. The accumulation of evidence from this and similar investigations 
may eventually lead to new typologies for dependence. Needless to say, much 
remains to be done before this goal can be achieved. 

APPENDIX: AUXILIARY RESULTS 
The following result justifies combination procedure 4.3.4. 

PROPOSITION A.l. Suppose that X = Y^k=i w kZl, where the Z k are in- 
dependent A^(0, 1) random variables, w\ > W2 > • • • and E(X) = J2kLi w k < 
oo. Then x~ l logP(X > x) — ► — l/{2w\) as x— >od. 

Proof. Fix a < l/(2w\). Since log(l — x) < x/(l — x), one has 
log{E(e QX )} = ~f] log(l - 2aw k ) < < oo. 
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Hence, by Markov's inequality one obtains P(X > x) < e - ax E(e aX ) and 
so limsup^oo x~ l logP(X > x) < —a. Letting a — > l/(2w\), one obtains 
limsupj-^oo x" 1 logP(X > x) < — l/(2w\). This concludes the proof since in 
view of large deviation results for Gaussian variables, one also has 

liminf x~ 1 logP(X > x) > liminf x~ l \ogP(w\Z 2 > x) = — — — . □ 

X— >oo x— too 2u)\ 

The following result is instrumental in establishing Proposition 6.1. 

Proposition A. 2. Let (w k ) be a positive sequence with J2fc=i w k < 00 
and (nk) be a sequence such that J2T=i w kf J 'k < oo. Furthermore, let (Z k ) be 
a sequence of independent jV(0, 1) random variables. For any 5 > 0, set 

oo 
k=l 

Then {P(X$ > x) - P(X > x)}/5 2 — > Y,V=i w kl^\ hk(x) as 5 — > 0, where h k 
is a density whose associated characteristic function /(t,0)/(l — 2iwkt) = 
(1 — 2iuikt)~ 1 rj^=i(l — 2iuijt)~ 1 / 2 is that of Xq -\-WkX2' ^ n which the sum- 
mands are taken to be independent. 

Proof. It follows from the definition of X$ that 



f(t,S) =E(e itXs ) = J] fj(w k t,6n k ) =t{t)e~ 252t2 ^ (t) e^M+^W 
fc=i 

where 



fc=i fe=i 

oo oo 



ft=l fc=l 

Note that and t 2 £{t) are integrable, «i is bounded, Ki(t)/t is bounded 
for i G {2, 3} and that as t — > 0, K2(t)/t — ► 1/36 and K%{t)/t — > J2T=i w kf J 'k- 
Next, from the Gil-Pelaez representation, one has 

1 r+oo 

P{X & >x)- P(X >x) = — t~ l Im{e~* tx f(t, 5) - e~* tx f(t, 0)} dt. 

ZTT J-oo 

Note that 

t- 1 lm{e- itx f(t, 6)} = t- 1 e(t)e- 2,52 ' 2KlW sin{ K2 (t) + S 2 K 3 (t) - tx}. 
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As a result, (5 2 t) 1 Im{e ltx f(t,5) — e ltx f(t,0)} can be decomposed as 
A^t^^^t) + A 2 (t,5)£(t), where 

Ai(i,tf) = (^^J-^e- 2 * 9 * 3 " 1 ^) - l}sin{K 2 (i) + S 2 K 3 (t) - tx} 

and 

A 2 (t,5) = ((5 2 t) _1 [sin{K 2 (0 + <5 2 K 3 (t) - tx} - sin{re 2 (i) - tx}}. 

Both terms are bounded and converge as 5 — > 0. Their limits are A\ (t, 0) = 
— 2t~ 1 Ki(t) sin{K2(t) — tx} and A 2 (t,0) = t~ x Ks(t) cos{n 2 (t) — tx}, respec- 
tively. An application of Lebesgue's dominated convergence theorem thus 
yields 

/oo poo 
t- 1 lm{e- itx f(t,8)-e- itx f(t,0)}dt= / ip(t,x)dt, 
-oo J—oo 

where ip(t, x) = £(t)[i _1 K3(i) cos{n 2 (t) — tx} — 2tK\{t) sm.{n 2 {t) — tx}]. 
It is easy to check that ip can also be expressed as 

oo 

i/>(t,x) = Y / w kV 2 k Re{e- ltx f(t,0)(l - 2itw k y 1 }, 

k=l 

where Re(z) stands for the real part of any complex number z. 
Since £ is integrable, so is f(t, 0)/(l — 2itwk) and hence 

1 C CO oo 2 /*oo 

— / ^{t,x)dt = Y J ^ hL / Re{e-' tte />,0)(l-2^ fc )- 1 }^ 

oo 

k=l 

where hj. is the density of Xq + w k x 2 , whose summands are taken to be 
independent. This completes the proof. □ 
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